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Abstract Based on morphological character comparisons and molecular phylogenetic reconstruction, we assessed the 
validity of the ratsnake subspecies Elaphe carinata deqenensis. Phylogenetic relationships inferred from fragments 
of two mitochondrial genes (ND2 and ND4) and one nuclear gene (c-mos) revealed a well-supported and minimally 
divergent sister relationship between putative E. c. deqenensis and E. c. carinata. Morphological comparisons indicated 
that all diagnostic characters proposed in the original description of E. c. degenensis are unable to separate this taxon 
from E. c. carinata. On the basis of these preliminary investigations, we provisionally suggest that E. c. deqenensis is 


not a valid subspecies, and should henceforth be treated as a synonym of E. c. carinata. 
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1. Introduction 


The King ratsnake, Elaphe carinata (Günther, 1864), 
belongs to the family Colubridae (Pyron et al., 2011), 
and ranges through Vietnam, Japan and China. In China, 
this species is widely distributed across a large region 
from western Yunnan (Gongshan County.) in the west to 
Shanghai in the east, and from northern Guangdong in 
the south to Hebei in the north (Zhao, 2006). E. carinata 
is a large snake with adults ranging from 1500 mm to 
2000 mm in total length (Zhao et al., 1998). This snake 
occupies a variety of habitats at elevations ranging 
between 100 and 2500 m above sea level (Zhao, 2006). 
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Elaphe carinata was originally described as Phyllophis 
carinata based on a single specimen from China (Günther, 
1864). In addition to the nominate form, another two 
subspecies Æ. c. deqenensis and E. c. yonaguniensis were 
recognized (Uetz, 2011). E. c. yonaguniensis is endemic 
to the Ryukyu Islands, Japan and Taiwan, China. The 
subspecies E. c. deqenensis was described from Deqin 
County, Yunnan, China based on several specimens (Yang 
and Su, 1984). The taxon Æ. c. deqenensis is restricted to 
a small area in Yunnan (Deqin and Weixi counties) and 
is distinguished from E. c. carinata by body coloration, 
pattern, morphometrics and scalation (Yang and Su, 
1984; Yang and Rao, 2008). The validity of this taxon 
is currently disputed, with the arrangement accepted by 
Utiger et al. (2002) and Uetz (2011), but not by Zhao et al. 
(1998) and Zhao (2006). Despite these disagreements, no 
morphological or molecular study has been conducted to 
examine the systematics of this taxon. 

Three specimens of E. carinata were collected from 
the type locality of E. c. deqenensis (Deqin County, 
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Yunnan) during fieldwork in 2010 (Figure 1). These 
specimens have allowed investigation of the status of this 
subspecies on the basis of morphological comparisons 
and molecular phylogenetic reconstruction. 


Figure 1 Sample locality of E. c. deqenensis (indicated by a red 
dot). 


2. Materials and Methods 


2.1 Morphological comparisons Thirty morphological 
characters related to scalation, and body dimensions were 
recorded for the three specimens (YBU 10004-10006; 
two adult males and one juvenile male) (Appendix). 
Measurements of body and tail length were taken with 
a ruler to the nearest 1 cm. Ventral scale counts were 
recorded following Dowling (1951), while the preparation 
and description of hemipenes followed Dowling and 
Savage (1960). The type specimens of E. c. deqenensis 
were unavailable for examination as they could not be 
found in the museum in which they were deposited, 
therefore morphological data for these specimens were 


Table 1 Specimens used in the molecular phylogenetic analysis. 


obtained from the original description (Yang and Su, 
1984) and from Yang and Rao (2008). Comparative data 
from specimens from other E. carinata populations were 
obtained from published work (Wu et al., 1985; Huang, 
1990; Fan et al., 1998; Zhao et al., 1998; Yang and Rao, 
2008). 


2.2 Molecular phylogenetic reconstruction Tissue 
samples from two of the three putative E. c. degenensis 
specimens were used for DNA sequence generation, 
and comparative sequences from seven other species 
of Elaphe (sensu stricto) (Burbrink and Lawson, 2007) 
were obtained from GenBank (Table 1). Based on 
previous molecular studies examining the relationships 
of ratsnakes, several closely-related representatives of 
the subfamily Colubrinae were chosen as outgroup taxa 
(Burbrink and Lawson, 2007) (Table 1). 

Genomic DNA was extracted from 85% ethanol- 
preserved tissues using standard proteinase K and phenol- 
chloroform protocols (Sambrook and Russell, 2002). Two 
mitochondrial protein-coding gene fragments NADH 
dehydrogenase subunit 4 (ND4), NADH dehydrogenase 
subunit 2 (ND2), and a fragment of one nuclear protein- 
coding gene (the oocyte maturation factor c-mos) 
were amplified following Arévalo et al. (1994) (ND4), 
Burbrink and Lawson (2007) (ND2) and Slowinski 
and Lawson (2002) (c-mos). These amplified loci were 
sequenced using an ABI 3730 Genetic Analyzer (Applied 
Biosystems) following the manufacturer’s protocols. 
Novel sequences generated in the present study are 
deposited in GenBank (Accession Nos. JN799413-— 
JN799418) (Table 1). 

DNA sequences were aligned using the program Mega 
4.0 with default parameters (Tamura et al., 2007; Kumar 


Taxon Voucher No. Origin GenBank accession No. (ND2/ND4/c-mos) 
Elaphe bimaculata DQ902210/DQ902283/DQ902062 
E. carinata carinata LSUMZ 37012 China DQ902211/DQ902284/DQ902063 


YBU 10004/GP 1271 
YBU 10005/GP 1272 


E. carinata deqenensis 
E. carinata deqgenensis 


E. climacophora CAS 163993 
E. dione LSUMZ 45799 
E. quadrivirgata 

E. quatuorlineata LSUMZ 40626 
E. schrenki 

Pantherophis vulpinus (O) CAS 184362 
Rhinechis scalaris (O) LSUMZ 37393 
Stilosoma extenuatum (O) LSUMZ 40124 
Zamenis situlus (O) CAS 175031 


JN799417/JN799413/JN799415 
JN799418/JN799414/JN799416 


Yunnan, China 
Yunnan, China 


Japan DQ902212/DQ902285/DQ902064 
Russia DQ902214/DQ902287/DQ902066 
Japan DQ902228/DQ902300/DQ902078 
Turkey AY487028/AY487067/AY 486955 
DQ902233/DQ902302/DQ902082 
Ohio, USA DQ902238/DQ902306/DQ902089 
Spain AY487029/AY 487068/AY 486956 
Florida, USA DQ902245/AF138776/DQ902093 
Unknown DQ902234/DQ902303/DQ902083 


YBU: Yibin University; CAS: California Academy of Sciences, San Francisco; LSUMZ: Louisiana State Museum of Natural Science; GP: 
Authors’ catalogue numbers. Outgroups are followed by (O), and sequences generated in this study are shown in bold. 
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et al., 2008). Protein-coding fragments were translated 
into amino acid sequences using Mega 4.0 to check for 
premature stop codons, which may indicate amplification 
of pseudogenes (Zhang and Hewitt, 1996). Average 
divergence estimates between species and subspecies 
were calculated from each locus separately as well as 
from all three loci combined Mega 4.0 (Tamura et al., 
2007; Kumar et al., 2008). 

Bayesian inference (BI) and maximum parsimony 
(MP) were used to reconstruct phylogenetic relationships. 
For the Bayesian analysis, sequence data were partitioned 
by locus, resulting in three partitions (c-mos, ND4 and 
ND2). The simplest best-fit model of sequence evolution 
(Posada and Crandall, 1998; Posada and Buckley, 2004) 
was inferred under the Akaike Information Criterion 
(AIC) using MrModeltest v. 2.2 (Nylander, 2004). The 
Bayesian analysis was performed using MrBayes v. 
3.1.2 (Huelsenbeck and Ronquist, 2001; Ronquist and 
Huelsenbeck, 2003), with three independent runs each 
consisting of four Markov chains (three heated chains 
and a single cold chain) starting from random trees, and 
applying the appropriate models of sequence evolution. 
Each analysis was run for a total of 5 x 10° generations 
and sampled every 1000 generations. Stationarity was 
confirmed by plotting the likelihood against generation 
in the program Tracer v. 1.4 (Drummond and Rambaut, 
2007). MP trees were generated from unweighted 
characters using PAUP* 4.0b10 (Swofford, 2003), 
applying a heuristic search with 1000 random sequence 
addition replicates and tree bisection-reconnection (TBR) 
branch swapping. Support values for clades were obtained 
by bootstrapping (BS; Felsenstein, 1985) using 1000 
bootstrap pseudoreplicates. 


3. Results 


3.1 Description of morphology The everted hemipenes 
are single and bulbous, each with a simple sulcus that 
extends to the top of the organ. The basal 1/5 of each 
hemipenis is smooth, followed by an approximately equal 
segment covered in papillae. The top 1/5 of the organ is 
nude and smooth, and the middle 2/5 is calyculate, with 
the edges of each calyce ornamented with small spines 
(Figure 2). 

The following descriptions are made based on two 
adult males collected from Deqin County, Yunnan (Figure 
1). These snakes are large, stout-bodied animals. The 
snout-vent lengths (SVL) of the two specimens examined 
were 1450 mm and 1080 mm, with respective tail lengths 
(TL) of 310 mm and 280 mm, giving SVL/TL ratios of 
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4.67 and 3.86. Other characteristics were: Supralabials 8, 
with 3 being anterior to the orbit, 2 entering the orbit, and 
3 posterior to the orbit; infralabials 10 or 11, with the first 
pair in contact with one another, the first to fifth touching 
the anterior chin shields and the fifth and sixth touching 
the posterior chin shields; preoculars 2, postoculars 2, 
temporals 6 or 7, loreal 1; both males having 23—23—17 
dorsal scale rows, strongly keeled with the exception of 
the outer row on each side; ventral scales 213 + 2 and 218 
+ 2, subcaudals 71 and 88 pairs, anal plate divided. The 
coloration of both adult males is similar. Scales on the 
top of the head are rich brown, paling on the sides. This 
color extends onto the nape but rapidly becomes darker 
and less rich, and most body scales are a dull brown. The 
interstitial skin is black, with numerous irregular cream 
cross-bars. This cream color is also present on the bases 
of dorsal scales that coincide with the bars. The bars 
fade posteriorly to scattered cream flecks, and disappear 
altogether towards the vent. The underside of the head 
is yellow-white. Ventrals are basically white with 
innumerable tiny black speckles, which tend to be dense 
posteriorly (Figure 3). 


Figure 2 Sulcate (A) and asulcate (B) views of a hemipenis of E. 
carinata (YBU 10005) from Deqin County, Yunnan, China. 


The juvenile specimen measures 350 mm SVL and 
95 mm TL, with a SVL/TL ratio of 3.68. Ventral scales 
are 212 and subcaudals 99. The dorsum is light tan with 
four longitudinal stripes, two dorsolateral and two lateral, 
which originate immediately behind the head. The lateral 
stripes terminate around the position of the cloaca, while 
the dorsolateral stripes extend to the tip of the tail. Dark 
brown flecks are scattered sparsely and irregularly on the 
dorsum, with density decreasing caudad, and disappearing 
altogether approximately 2/3 of the way down the body. 
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Ventrals are consistent with the dorsal scales in coloration, 
but scattered black dots are present on both sides of the 
ventrals (Figure 4). 


Figure 3 Dorsal view of an adult male Æ. carinata (YBU 10005) 
from Deqin County, Yunnan, China. 


Figure 4 Dorsal view of the juvenile E. carinata (YBU 10006) 
from Deqin County, Yunnan, China. 


3.1 Molecular phylogeny The final molecular dataset 
consisted of 2292 DNA base pairs (bp): 1032 bp ND2, 
693 bp ND4, and 567 bp c-mos. 701 bp (30.6%) and 
614 bp (26.8%) were variable, and 486 bp (21.2%) and 
356 bp (15.5%) were parsimony-informative when the 
outgroups were included and excluded. We are confident 
that pseudogenes were not amplified as no insertions, 
deletions, or stop codons were detected in the protein- 
coding genes (Zhang and Hewitt, 1996). Uncorrected 
sequence divergence from interspecific comparisons 
(excluding outgroups) ranged between 12.4%-17.9% for 
ND2, 12.1%-17.3% for ND4, and 9.1%-—12.2% for the 
combined dataset. Intraspecific sequence divergences for 
E. carinata were much lower, ranging between 1.9%- 
2.8% for ND2, ND4 and the combined dataset (Table 2). 
The optimal models of sequence evolution identified 


by MrModeltest were as follows: GTR + I + T for ND4, 
HKY+ I +T for ND2 and HKY for c-mos. These models 
were implemented in our analyses, and the resultant BI 
50% majority-rule consensus tree is shown in Figure 5. 
All seven representatives of Elaphe were monophyletic, 
and the two individuals representing E. c. deqenensis 
were included with E. c. carinata in a well-supported 
clade (100% PP: Figure 5). 


Elaphe climacophora 


100/99 | 100/89 Elaphe bimaculata 


Elaphe dione 


Elaphe quatuorlineata 
80/- 
Elaphe carinata 


YBU 10004 


95 /— 
100/100 yey 10005 


Elaphe quadrivirgata 


Elaphe schrenckii 


Outgroups 


0.1 


Figure 5 Bayesian 50% majority-rule consensus tree inferred from 
the combined ND2, ND4 and c-mos data. Node support is provided 
in the format: Bayesian posterior probabilities / MP bootstrap 
values. 


The parsimony analysis produced a single tree with a 
length = 504, CI = 0.58, RI = 0.51, and RC = 0.30. The 
MP tree showed consistently high support values with BI 
tree in the clade including all putative species of Elaphe 
(99% BS) and the clade including the three specimens 
of E. carinata (100% BS) (Figure 5). The disagreements 
were indicated in several nodes which are very low in 
support indices in MP tree. Bootstrap values are reported 
on the Bayesian 50% majority-rule consensus tree in 
Figure 5. 


4. Discussion 


Traditionally, subspecies have been identified on the basis 
of morphological characters. However, morphological 
variation may result from phenotypic plasticity and 
adaptation to local environments, which may mislead our 
understanding of evolutionary relationships (Herrmann 
et al., 2004). DNA sequence data are, thus, commonly 
used for phylogenetic reconstruction, as is the case in this 
study. Both methods of phylogenetic tree reconstruction 
showed that the specimens from Deqin currently 
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Table 2 Average interspecific and intraspecific divergence estimates (% Tamura-Nei distances, with gamma correction). Above diagonal: 
Distances calculated from ND2/ND4; Below diagonal: Distances calculated from the combined data set (ND2, ND4 and c-mos). 
Comparisons between E. c. deqenensis and E. c. carinata are shown in bold. The distances calculated from c-mos are not included due to 


their low variation. 


c. deqenensis bimaculata c. carinata climacophora dione quadrivirgata quatuorlineata schrenckii 
c. deqenensis 13.9/15.1 2.8/1.9 17.2/14.9 15.4/17.2 16.7/14.1 14.8/15.7 15.0/14.9 
bimaculata 10.8 14.5/14.1 16.9/14.7 12.4/12.4 15.1/15.7 12.9/13.7 14.5/16.2 
c. carinata 2.0 10.7 17.5/13.6 15.0/16.7 18.0/12.1 16.3/14.2 16.4/14.2 
climacophora 12.1 11.7 11:7 17.9/15.1 16.5/14.0 15.1/14.2 17.0/15.0 
dione 12.1 9.1 11.7 12.2 15.2/15.0 14.6/15.4 14.9/17.3 
quadrivirgata 11.5 11.3 11.4 11.4 11.2 14.9/15.7 13.4/14.9 
quatuorlineata 11.2 9.7 11.3 10.8 10.9 Lid 13.5/16.3 
schrenckii 11.1 11.2 11.4 11.8 113 10.3 10.6 


recognized as E. c. degenensis form a well-supported 
clade with Æ. c. carinata, and this clade is nested within 
the genus Elaphe (Figure 5). Based on our small sample, 
the genetic divergence between the two subspecies of 
Elaphe carinata is approximately one quarter that between 
the Elaphe species included in this study (Table 2). When 
compared with other studies of colubroid snakes, the 
genetic divergence within E. carianta is low (Guicking 
et al., 2008; Guo et al., 2011). For example, for ND4, the 
intraspecific genetic divergence of Protobothrops jerdonii 
varies between 2.8%—4.7% (Guo et al., 2011), which is 
much higher than that of E. carinata; but the interspecific 
divergence of Protobothrops ranges between 3.5%-13% 
(Liu et al., 2012), which is consistent with those of 
Elaphe. 

The hemipenes of E. c. degenensis and E. c. carinata 
appear to be morphologically similar based on the 
descriptions of the specimens from Fujian, China by 
Pope (1935) and the specimens from Hunan, China by 
Zhang et al. (1984). However, morphological data from 
our putative E. c. deqenensis specimens do not fit with 
the meristic and morphometric characters proposed for 
subspecies delineation in E. carinata by Yang and Su 
(1984), and Yang and Rao (2008). A detailed comparison 
of our results with those of Yang and Su (1984), and Yang 
and Rao (2008) is presented below, where the phrases in 


Table 3 Body measurements of E. carinata from China. 


bold are direct quotations from these two papers. 

Lower ratio of SVL to TL: After the original description 
(Yang and Su, 1984), Yang and Rao (2008) added the ratio 
of SVL to TL as a diagnostic character for E. c. deqenensis. 
Yang and Rao (2008) proposed that in degenensis this ratio 
is lower (3.1 in males and 2.8 in a single female) than in £. 
c. carinata (4.36 in males and 4.84 in females). However, 
the ratios for the two male and one juvenile degenensis 
examined in this study (ranging from 3.68 to 4.67) were 
higher than those reported by Yang and Su (1984), and 
Yang and Rao (2008). In addition, some specimens 
belonging to E. c. carinata have lower ratios of SVL to 
TL. For example, a male from Guizhou has a ratio 3.17 
(Wu et al., 1985), and two females from Shanxi have ratios 
of 2.92 and 2.84 (Fan et al., 1998) (Table 3). Though our 
sampling is not sufficient to investigate the distribution 
of this character in the two putative subspecies, it is clear 
that it cannot serve as a diagnostic character for their 
separation. 

Ventral scales are less than 214: Yang and Su (1984), 
and Yang and Rao (2008) proposed that deqenensis has 
fewer than 214 ventral scales (210.6 in males and 214 in 
a female), while E. c. carinata has more than 214. Based 
on published data it appears that putative E. c. carinata 
specimens from Zhejiang, Sichuan, Guizhou, Shanxi, 
Gansu and Guangxi often exhibit ventral scale counts 


Locality Sex SVL (mm) TL (mm) SVL/TL References 
Deqin, Yunnan 288 850, 880 280, 320 2.75, 3.04 Yang and Rao (2008) 
Deqin, Yunnan 288 1450, 1080 310, 280 4.67, 3.86 This study 
Gongshan, Yunnan 288 1070, 1106 240, 260 4.25, 4.46 Yang and Rao (2008) 
Songtao, Guizhou 12 88 940-1560 238—436 3.17-4.67 Wu et al. (1985) 
Zhejiang 1033 7 90-1430 195-356 3.83-4.50 Huang et al. (1990) 
Tengchong, Yunnan 799 1112-1540 229-327 4.64-5.07 Yang and Rao (2008) 
Songtao, Guizhou 16 99 934-1666 220-310 3.88-6.26 Wu et al. (1985) 
Zhejiang 822 820-1535 200-360 4.10-5.24 Huang et al. (1990) 
Shanxi 299 1022, 1010 350, 356 2.92, 2.84 Fan et al. (1998) 
Deqin, Yunnan J 350 95 3.68 This study 
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of less than 214 (Table 4). We also found that one of our 
putative degenensis specimens (YBU 10005) had 218 
ventral scales. Additionally, the number of subcaudals 
in degenensis varies between 71 and 94 (Yang and Rao, 
2008; this study), which is within the range of variation 
in E. c. carinata (69—99, though one specimen from 
Zhejiang had a count of 57; Table 4). 

Body scales at neck are 25: Yang and Su (1984) and 
Yang and Rao (2008) proposed that degenensis has 25 
rows of body scales on the neck, and they regarded this 
as a diagnostic character for degenensis. However, this 
character is considerably variable in most snakes (e. g., 
Protobothrops jerdonii, personal observation), thus it is 
seldom used for taxon delimitation as it is of little use 
in discriminating between snake species and subspecies 
(Zhao et al., 1998). 

Different body coloration and pattern: The coloration 
and patterning of our putative degenensis specimens were 
consistent with the descriptions provided by Yang and Su 
(1984), and Yang and Rao (2008), in which degenensis 
lacks yellow coloration on the body, has narrow black 
cross-bands at the anterior of the body, and no reticulate 
pattern posteriorly or on the tail. E. carinata displays 
intraspecific color and pattern polymorphism that appears 
to be associated with habitat (Zhao et al., 1998; Zhao, 
2006; personal observation) and with age. Intraspecific 
variation in color and pattern in snakes can be extreme, 
for example, in Pantherophis obsoletus (Burbrink et al., 
2000), Himalayophis tibetanus (Gumprecht et al., 2004; 
Vogel, 2006), Protobothrops jerdonii (Gumprecht et al., 


2004; Vogel, 2006; Guo et al., 2009), Opisthotropis cheni 
(Li et al., 2010). Such variation has often formed the basis 
for description of species and subspecies (e. g., Orlov 
and Helfenberger, 1997, but see Burbrink et al., 2000). 
However, comprehensive evolutionary investigations 
(particularly those using molecular phylogenies) often 
do not support taxonomic arrangements based solely on 
coloration and pattern (Burbrink et al., 2000; Tillack et 
al., 2003; Guo et al., 2009). For example, Pantherophis 
obsoletus was historically regarded as comprising up to 
eight subspecies defined primarily on the basis of adult 
color and pattern. However, molecular phylogenies 
coupled with detailed morphological analyses 
demonstrated that these subspecies did not represent 
distinct evolutionary lineages (Burbrink et al., 2000). Thus, 
Burbrink et al. (2000) proposed that it was dangerous to 
recognize subspecies based on few characters, especially 
those associated with color and pattern. 

Although the molecular results were not sufficient 
to conclude whether the subspecies E. c. degenensis 
was a valid evolutionary unit or not, the morphological 
comparison indicated that the putatively diagnostic 
morphological characters could not be used to separate E. 
c. deqenensis from E. c. carinata. Thus, we provisionally 
suggest that the populations of E. carinata from 
Deqin and Weixi counties of China do not represent 
distinct taxonomic entities, and should be regarded 
as members of the typical race (E. c. carinata). More 
extensive geographic sampling and sequences from 
additional unlinked genetic markers will be necessary to 


Table 4 Ventral (Vs) and subcaudal (Sc) scale counts for E. carinata from China. 


Locality Sex Vs Sc References 
Deqin, Yunnan 10 3¢ 206-213 72-94 Yang and Rao (2008) 
Deqin, Yunnan 28g 213,218 71, 88 This study 
Gongshan, Yunnan 243 217-221 83 Yang and Rao (2008) 
Guizhou 5243 194-228 81-103 Wu et al. (1985) 
Zhejiang 1033 213-222 85-99 Huang et al. (1990) 
Guangxi 13 203 97 Zhao et al. (1998) 
Zhejiang 23 88 212-227 75-102 Zhao et al. (1998) 
Sichuan 38g 208-218 87-93 Zhao et al. (1998) 
Shaanxi 18 214 73 Zhao et al. (1998) 
Deqin, Yunnan 19 214 92 Yang and Rao (2008) 
Deqin, Yunnan J 212 92 This study 
Tengchong, Yunnan 799 218-224 74-83 Yang and Rao (2008) 
Guizhou 60 22 186-226 73-95 Wu et al. (1985) 
Zhejiang 899 216-220 57-92 Huang et al. (1990) 
Zhejiang 15 99 198-219 73-89 Zhao et al. (1998) 
Shanxi 299 208-209 84-86 Fan et al. (1998) 
Guangxi 19 219 86 Zhao et al. (1998) 
Yizhang, Hunan 1:9 220 69 Zhao et al. (1998) 
Sichuan 299 213-214 82-87 Zhao et al. (1998) 
Gansu 299 211 80-81 Zhao et al. (1998) 
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comprehensively investigate the possible existence of 
cryptic taxa within the polymorphic E. carinata. 
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Appendix 


Morphological characters used in the morphological analysis, and 
their abbreviations are provided below: 


1. Scalation 

BSCK: Keeling of body scales (mid body), recorded as 0: none; 0.5: 
weak; 1: strong. 

KHEADSC: Keeling of the scales on the back of the head. 

KTEMP: Keeling of the temporal scales. 

SC: Number of pairs of subcaudal scales (any unpaired scales are 
treated as a pair). 

SUBLAB: Average number of sublabials on the left and right hand 
side. 

SUPLAB: Average number of supralabials on the left and right hand 
side. 

VS: Number of ventral scales (VS), not including anal scale, 
recorded by the Dowling (1951) method (i. e., the first VS is the one 
which contacts the first dorsal scale row on both sides). 

PREOC: Number of preocular scales. 

POSTOC: Number of postocular scales. 

Anal: Divided or fused. 

NSC: Numbers of the single subcaudal scale. 

Loreal: Numbers of the loreal. 


2. Scale reduction formula 

These are recorded as a series of characters, each referring to a 
specific reduction. Each position will have two characters, the 
dorso-ventral (DV) position of the reduction (the lowest of the two 
merging scale rows), and the ventral scale (VS) position (counted 
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from the head), which is the ventral scale to which the scale 
reduction traces diagonally. Before analysis, the VS position was 
transformed into the percentage of the total number of ventral scales 
(%VS), to control for variation. (Sbucaudal scale position: SC). 
VS21tol9: Ventral scale position of the reduction from 21 to 19 
scale rows. 

VS19tol7: Ventral scale position of the reduction from 29 to 27 
scale rows; and then SC6to4, SC8to6, DV21tol9, DV19tol7, 
DV8to6, DV6to4. 


3. Body dimensions 

All measurements are made on the right side of the head only unless 
it was damaged, in which case they were done on the left. 

DEYE: Diameter of the eye measured between the edges of the 
surrounding scales. 

EYE2NOS: Distance between the eye and the nostril, measured 
between the suture between the second and third preocular (from 
the bottom) and the inner edge of the nostril. 

LHEAD: Length of the head, between the tip of the snout to the 
posterior edge of the lower jawbone. 

LSUPOC: Length of the supraoculars measured. 

WSUPOC: Width of the supraoculars measured. 

SVL: Distance between the tip of the snout and the cloaca. 

TAIL: Distance between the anterior edge of the first subcaudal 
scale and the tail tip. 

WMT: Width of the top rostal. 

WMB: Width of the bottom rostal. 

WHEAD: Width of the head measured between the outer edges of 
the supraoculars. 

WSUPOC: Width of the supraoculars measured at the widest part. 


